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Abstract. We report on an investigation into stellar evolution and nucleosynthesis in the low and extremely low metallicity 
regime, including models of stars with a pure Big Bang composition (i.e. Z = 0). The metallicity range of the extremely metal- 
poor (EMP) models we have calculated is —6.5 < [Fe/H] < —3.0, whilst our models are in the mass range 0.85 < M < 3.0M,t). 
Many of the EMP and Z = models experience violent evolutionary episodes not seen at higher metallicities. We refer to these 
events as 'Dual Flashes' since they are characterised by peaks in the hydrogen and helium burning luminosities occurring at 
roughly the same time. Some of the material processed by these events is later dredged up by the convective envelope, causing 
significant surface pollution. These events have been reported by previous studies, so our results reaffirm their occurrence 
- at least in ID stellar models. The novelty of this study is that we have calculated the entire evolution of the Z = and 
EMP models, from the ZAMS to the end of the TPAGB, including detailed nucleosynthesis. We have also calculated the 
nucleosynthetic yields, which will soon be available in electronic format. Although subject to many uncertainties these are, as 
far as we are aware, the only yields available in this mass and metallicity range. In this paper we briefly describe some of the 
results in the context of abundance observations of EMP halo stars. This work formed part of SWC's PhD thesis (completed 
in March 2007) and a series of subsequent papers will describe the results of the study in more detail. 
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MOTIVATION 

The discovery of extremely metal-poor stars (EMPs) in 
the Galactic Halo has renewed interest in the theoreti- 
cal modelling of Population III and low-metallicity stars. 
Most of these stars show abundances that conform to a 
simple Galactic chemical evolution line (see Figure [TJ. 
However a subset of the EMP stars have been observed 
to contain large amounts of carbon. These C-rich EMPs 
(CEMPs) make up a large proportion of the EMP popu- 
lation (~ 10 — -> 20%; see eg. fljj] ) . This population is also 
highlighted in FigureQ] Apart from carbon the EMP stars 
also display variation in a range of other elements (see 
0]] for a review of the observations). A number of theo- 
ries have been proposed to explain the various patterns, 
ranging from pre-formation pollution via Pop III super- 
novae (eg. J2l, Hi) to self-pollution through peculiar evo- 
lutionary events (eg. Jj], ||5|]) to binary mass transfer (eg. 

In the current study we have undertaken a broad ex- 
ploration of EMP stellar evolution and nucleosynthe- 
sis in the low and intermediate mass regime. Our study 
expands on the previous work in the field. In particu- 
lar it includes (for the first time) full evolutionary and 
nucleosynthesis calculations from ZAMS to the end of 
the thermally-pulsing AGB phase (TPAGB), as well as 
chemical yields for 74 nuclear species. With this homo- 



geneous set of models we hope to shed some light on 
whether or not ID stellar models can explain some of the 
EMP halo star observations. 



METHOD 

Our simulations were performed utilising two numerical 
codes - a stellar structure code and post-process nucle- 
osynthesis code. 

The stellar structure code used was the Monash/Mount 
Stromlo STAR code (MONSTAR see eg. 01, HI). 
We do not describe the code in detail here but briefly 
note a few key points. The code is largely a standard 
ID code that utilises the Henyey-matrix method (a mod- 
ified Newton-Raphson method) for solving the stellar 
structure equations. Opacities have been updated to those 
from 1 16] (for mid-range temperatures) and 11711 (for 
low temperatures). For the present study the instanta- 
neous convective mixing routine was replaced by a time- 
dependent (diffusive) mixing routine (similar to that de- 
scribed by [18]). This change was necessary due to the 
violent evolutionary events that occur in models of EMP 
stars. Convective boundaries were always defined by the 
Schwarzschild criterion - ie. no overshoot was applied. 
Thus the extension of convective zones in all the models 
are conservative. A key problem with modelling EMP 
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FIGURE 1. Carbon abundance observations versus [Fe/H] 
for metal-poor stars compiled from the literature. The dotted 
line marks [C/Fe] = 0.0. The normal EMPs lie on this Galactic 
chemical evolution line, whilst the CEMPs (which we define 
as having [C/Fe] > +0.7) lie above this. The two most Fe-poor 
stars known, HE 0107-5240 ([Fe/H] = -5.3, [7]) and HE 1327- 
2326 ([Fe/H] = -5.4, H]) are labelled. We note that this sample 
is biased towards CEMPs. Observational data sets are from [9], 

oa, o, m, and m. 



stars is the unknown driver(s) of mass loss. The domi- 
nant theory is that mass is lost from red giant envelopes 
through radiation pressure acting on grains. Thus, in the 
EMP or Z = regime, mass loss is thought to be neg- 
ligible. The MONSTAR code uses empirical mass loss 
formulae (the formula from [19] during the RGB and 
that of [20] during the AGB). For this study we chose to 
retain the standard treatment, for the following reasons. 
The first is that not much mass is lost during the RGB 
in these EMP models. This is because their RGB phases 
are much shorter (or non-existent), since He is ignited 
much earlier than in stars of comparable mass at higher 
metallicities. Indeed, the luminosity at the tip of the RGB 
can be up to 1 order of magnitude lower in EMP models 
(see eg. [21]). Thus the RGB mass loss is largely negligi- 
ble, and it is the AGB mass loss that needs to be handled 
properly. As described in the next Section these models 
all experience some sort of polluting episode - and al- 
ways before the AGB phase. This has the consequence 
that the surface of the AGB models usually have metal- 
licities approaching that of the LMC or even Solar (as 
defined by Z rather than Fe - they are still metal poor in 
terms of Fe). Thus, since the stellar surfaces have (some 
of) the ingredients needed to form grains, we argue that 
using the standard mass loss formula given by 112011 is 
warranted. We also note that metallicity is also indirectly 
taken into account by the mass loss formulae, since they 
depend on bulk stellar properties (such as radius, lumi- 
nosity, pulsation period), which vary significantly with 
metallicity. 



The nucleosynthesis calculations were made with 
the Monash Stellar Nucleosynthesis code (MONSN). 
MONSN is a post-process code. As input it takes the 
key structural properties of each hydrostatic model from 
the MONSTAR code (eg. density, temperature, convec- 
tive velocities). It solves a network of 506 nuclear re- 
actions involving 74 nuclear species (see eg. rf22ll . 1I23I1 . 
l24ll for more details on this code). Our grid of models 
covered the mass range: M = 0.85, 1.0, 2.0,3. 0M S and 
the metallicity range: [Fe/H] = -6.5, -5.45, -4.0, -3.0, 
plus Z = 0. 



THE DUAL FLASH EVENTS 

It has long been known that theoretical models of Z = 
stars (and EMP stars) undergo violent evolutionary 
episodes not seen at higher metallicities. This was first 
suggested by B21I1 and confirmed by later calculations 
(see eg. J25[|, 112611 . rf27ll ). These episodes occur at dif- 
ferent evolutionary stages in stars of different mass and 
metallicity (see eg. |4fl). The most severe of these evo- 
lutionary events occurs during the core He flash of low- 
mass stars (with [Fe/H] < —2.5, see J3l). In this event 
the normal flash-driven convective zone breaks out of 
the He-rich core. Thus H-rich material is mixed down 
into regions of high temperature. Processed material is 
also mixed upwards. The proton-rich material burns at 
a very high rate, amounting to a secondary flash - a H- 
flash. This flash reaches luminosities comparable to the 
core He flash itself and occurs within the same time- 
frame as the He flash. Thus we refer to the combination 
of these events as a 'Dual Core Flash' (DCF). We note 
that this event has also been referred to as Helium Flash 
Induced Mixing (HEFM, J2I]) and Helium Flash-Driven 
Deep Mixing (He-FDDM, |6|]). A similar event occurs in 
stellar models of higher mass and higher (although still 
very low) metallicities. In these cases it is the normal 
AGB shell He flash-driven convection zone that breaches 
the H-He discontinuity. This occurs during the first few 
pulses of the TPAGB phase. Again a H-fiash is induced 
during the evolution of the He-flash, so we refer to this 
event as a 'Dual Shell Flash' (DSF). Both the DCF and 
DSF events have consequences for the surface composi- 
tion of the star since, in both cases, the convective enve- 
lope subsequently deepens and mixes up the (processed) 
material overlying the H-burning shell. Presently there is 
reasonable consensus that these events do indeed occur 
(at least in stellar models!), although we note that not ev- 
ery study has found them to (eg. B29I1 ). In Figure [2] we 
display an example calculation of a DCF from one of our 
models. 
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FIGURE 2. An example of one of the Dual Core Flashes 
that occurred in some of our low metallicity models. In this 
case it is the 1 M e model with [Fe/H]= —6.5. Convective 
regions are shown by grey shading whilst the mass location 
of the edge of the H-exhausted core is shown by the solid line 
(blue). The large inward movement of the H-He boundary (H- 
Ex. core), due to the ingestion of protons into the HeCZ, can 
be seen (AM f ). The beginning of the post-DCF dredge-up is 
indicated. This dredge up (mainly of CNO isotopes) increases 
the surface metallicity to super-solar values (although [Fe/H] 
remains unchanged). 

CATEGORISATION OF YIELDS 

In Figure [3] we show the results from a nucleosynthesis 
calculation for one of our models. In this model it is the 
Third Dredge-Up (3DU - the periodic dredging up of He 
burning products into the convective envelope) and Hot 
Bottom Burning (HBB - hydrogen burning at the bottom 
of the convective AGB envelope) that define the yield of 
the star. Indeed, the chemical signature arising from the 
DSF occurring at the start of the AGB is totally erased 
by these normal AGB evolutionary episodes. However 
this is not always the case. In Figure [4] we summarise 
the pollution episodes over the whole grid of models. We 
group the yields into four categories, defined by the evo- 
lutionary events/phases that dominate the chemical sig- 
nature in the yields: Group 1 yields are dominated by 
the DCF events, Group 2 are dominated by DSF events, 
Group 3 are dominated by 3DU+HBB, whilst Group 4 
show no surface pollution during their hydrostatic evo- 
lution (but may explode as Type 1.5 supernovae). Two 
key features visible in this figure are (1) members of the 
DCF group have polluted surfaces during the horizontal 
branch phase (HB, the core He burning phase) and (2) 
both the DCF and DSF groups, which are of low mass, 
have polluted surfaces during the AGB - despite the lack 
of 3DU. Thus our models predict a greater proportion of 
C-rich stars at extremely low metallicity, as these events 
do not occur at higher metallicities. 
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FIGURE 3. The time evolution of the surface composition 
in the [Fe/H] = —5.45, 3 M model, for selected species. 
This model fits into our Group 3. The rich nucleosynthesis 
arising from 3DU and HBB is seen in the right-hand panel. In 
particular the CN cycling product 14 N dominates during most 
of the AGB. The I2 C/ 13 C and C/N ratios quickly approach 
equilibrium values once the (strong) HBB starts. This chemical 
signature by far dominates that of the DSF in this case. 

COMPARISONS WITH OBSERVATIONS 

In Figure [5] we compare the carbon yields from our 
entire grid of models with the observed [C/Fe] abun- 
dances in EMP halo stars. It can be seen that the yields 
from our models are all C-rich. Moreover, the yields at 
[Fe/H] = —3, for which there are the most observations, 
agree well with the observations. At [Fe/H] = —4 there is 
a fair agreement with the observations, although there are 
less observations to compare with. At [Fe/H] = —5.45 
the yields show somewhat more C than is observed, al- 
though there are only two stars observed at this metal- 
licity. An interesting feature of this diagram is that the 
the model yields predict [C/Fe] to continue increasing 
towards lower metallicities. Furthermore, taking into ac- 
count the evolutionary stage at which the surface pollu- 
tion is gained in the lower mass models (M = 0.85 and 
1.0 Mp) - ie. the Dual Core Flash events rather than 
the AGB - the models also predict a higher proportion 
of C-rich stars at lower and lower metallicities. This is 
due to the fact that these stars already have polluted sur- 
faces during the HB stage - which has a lifetime roughly 
1 order of magnitude longer than the AGB phase. The 
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FIGURE 4. Displayed in each HR diagram is a representa- 
tive example from our grid of models for each pollution group. 
Grey (red) lines indicate phases of the evolution in which the 
surface is strongly polluted with CNO nuclides (from the DCF, 
DSF or 3DU events). Black (blue) lines indicate that the surface 
still retains the initial metal-poor composition. Evolutionary 
stages and pollution sources are marked, as are the mass and 
metallicity ranges of each group. Question marks indicate un- 
known upper boundaries (due to the limited mass range of the 
current study). 



CEMPs also show interesting behaviour in other elemen- 
tal abundances. In Figure [6] we compare the integrated 
yields from our [Fe/H] = —5.45 models with the ob- 
served abundance patterns of the two most metal-poor 
halo stars. Apart from N these two stars show similar 
abundance patterns. They both have roughly the same C 
and O abundances and, although offset to lower abun- 
dances, the Na and Mg abundances follow the same pat- 
tern. Nitrogen is 2 dex more abundant in HE 1327-2326. 
The dominance of N over (enhanced) C and O in HE 
1327-2326 is reminiscent of the pattern that CNO cy- 
cling plus 3DU creates. This pattern is also seen in the 
yields of the 2 M and 3 M models. In these models 
it is 3DU+HBB that dominates the chemical signature 
of the yields. It may be possible that HE 1327-2326 has 
been polluted by an intermediate mass AGB star that had 
undergone 3DU and partial (or lower temperature) HBB. 
Na and Mg are very much overproduced in these models 
compared to the observations but, again, lower tempera- 
ture or incomplete HBB may account for this. HE 0107- 
5240 on the other hand is better 'matched' (it is still far 
from ideal) by the 1 M© yield. The yield of this model is 
dominated by the DSF event at the start of the AGB. In 
this case N and Na are overproduced in the model rela- 
tive to the observations. 
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FIGURE 5. Comparing the yields from all our models with 
observations of EMP stars. See Figure Q] for the observational 
data sources. Here we have colour-coded the observations 
into [C/Fe]-rich (black dots) and [C/Fe] -normal (grey/magenta 
dots), where [C/Fe]-rich is defined by [C/Fe] > +0.7. The short 
horizontal lines indicate the starting composition of the models 
(in this case they are all at [C/Fe] ~ 0, except for the Z = mod- 
els). We have plotted the Z = model yields at [Fe/H] — —8 
for comparison. The two most metal-poor stars can be seen at 
[Fe/H] ~ —5.5 (they are both C-rich). An upper envelope to the 
pollution of the models - and the observations - is marked by 
the dotted line at [C/H] — +1.0. The yields from our models 
are colour- and shape-coded to highlight the different episodes 
that produced the bulk of the pollution in each yield. Numbers 
beside each yield marker indicate initial stellar mass, in M@. 
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FIGURE 6. Comparing the abundance pattern in our 
[Fe/H] = —5.45 yields (open circles, initial stellar masses as 
indicated) with observations of stars with the same metallicity 
(closed squares). The abundance determinations in these stars 
are from @] (HE 1327-2326) and Q] (HE 0107-5240). 



More comparisons and analysis will be reported in 
future papers arising from this thesis. 



SUMMARY 

In agreement with most previous studies we find that 
many of the EMP and Z = models experience violent 
evolutionary episodes not seen at higher metallicities. We 
refer to these events as 'Dual Flashes' (DFs). 

The models predict an increased number of carbon- 
rich stars at the lowest metallicities. This is mainly due 
to the extra pollution provided by the DF events. This 
concurs well with the observations, which show the pro- 
portion of CEMP stars in the Galactic halo to be higher 
at lower metallicities. Although not discussed here we 
also note that we found the pollution arising from the 
DF events to be simultaneously C- and N-rich, as also 
observed in the CEMP stars. This contrasts with the pol- 
lution expected from 3DU at low mass, which would be 
(mainly) C-rich. Furthermore, the models predict that the 
proportion of CEMP stars should continue to increase at 
lower metallicities, based on the results that the low mass 
EMP models already have polluted surfaces by the HB 
phase, and that there are more C-producing evolution- 
ary episodes at these metallicities. We also compared the 
chemical pollution arising from the models with the de- 
tailed abundance patterns available for some of the most 
metal-poor CEMP stars, and found mixed results. 

We note that all these calculations contain many uncer- 
tainties. These include the unknown mass-loss rates, un- 
certain nuclear reaction rates, and the treatment of con- 
vection. In the case of the DF events a further uncertainty 
is the possibility that full fluid dynamics calculations are 
really needed to model these violent episodes. 

The models and yields from this thesis will be de- 
scribed in more detail in a series of future papers. For 
interested readers a full version of the thesis is available 
on SWC's webpage. 
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